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1 Introduction 



The notion of distribution amplitudes (DAs) [T] refers to matrix elements of nonlocal light-ray 
operators sandwiched between the hadron state and the vacuum. The physical interpreta- 
tion of DAs is transparent in the infinite momentum frame [2], in which case DAs correspond 
to momentum-fraction distributions of partons in a hadron at small transverse separation. 
Equivalently, DAs can be related to the transverse momentum integrals of the hadron's 
Bethe-Salpeter wave functions that appear e.g. using the formalism of the light-cone quan- 
tization Schematically, 



The natural application of DAs in phenomenology are exclusive hard processes with 
large momentum transfer which can be calculated using factorisation methods. Apart from 
meson/baryon electromagnetic form factors, this includes a large class of phenomenologically 
very interesting B meson decays, such as weak decay form factors [1], the non-leptonic decays 
B —>■ M1M2, where Mi is a light meson [5], and also rare radiative and semileptonic decays, 
B M7 [6] and B —>■ Mi'^i~ [7], which involve flavour- changing neutral currents and are 
heavily suppressed in the Standard Model, but sensitive to new-physics effects. For these 
processes, the shape of the DAs is very important: in B ^ p'j vs. B — ^ K*'y, for example, 
the size of SU(3) breaking in the relevant DAs is presently the dominant source of theoretical 
uncertainty [8]. All these decays will be studied in detail at the forthcoming LHC, which 
makes the detailed investigation of the various leading and higher-twist DAs both timely and 
relevant. 

The crucial point and main technical difficulty in the construction of higher-twist DAs is 
the necessity to satisfy the exact equations of motion (EOM), which yield relations between 
physical effects of different origin: for example, using EOM, the contributions of orbital 
angular momentum in the valence component of the wave function can be expressed (for 
mesons) in terms of contributions of higher Fock states. An appropriate framework for 
implementing these constraints was developed in Ref. [S]: it is based on the derivation of 
EOM relations for non-local light-ray operators [10], which are solved order by order in the 
conformal expansion; see Ref. [11] for a review and further references. In this way it is 
possible to construct self-consistent approximations for the DAs, which involve a minimum 
number of hadronic parameters. Another approach, based on the study of renormalons, was 
suggested for twist 4 in Refs. [HI [13]: this technique is appealing as it allows one to obtain an 
estimate of high-order contributions to the conformal expansion which are usually omitted. 
In Ref. IH] we have generalized this approach to include SU(3)-breaking corrections and 
have shown how to combine renormalon-based estimates of "genuine" twist-4 effects with 
meson mass corrections. 

In a series of previous papers, Refs. [151 [IB] , we have developed the corresponding for- 
malism for vector mesons. In Ref. [T7] the analysis of twist-3 DAs of light vector mesons, p, 
K* and 0, was completed, including all SU(3) and G-parity-breaking effects. In this paper. 
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we put together the last missing pieces, completing the study of the corresponding DAs of 
twist 4, with the main emphasis on the calculation of SU(3)-breaking effects in the relevant 
hadronic matrix elements. These corrections come from different sources: 

• SU(3) breaking of hadronic parameters: these effects are known for twist-2 and -3 
parameters, see Refs. [T7l [TS] , but have not been studied for twist-4 parameters before; 

• G-parity- breaking parameters: these are of parametric order — and vanish in the 
limit of equal quark mass, i.e. for p and 0. For twist-2 DAs, they have been calculated, 
to lowest order in the conformal expansion in Refs. [T8l [T9l EHl El], and for twist-3 
DAs in Ref. they are unknown for twist-4 DAs; 

• explicit quark-mass corrections in m^img to DAs and evolution equations: these affect 
only higher-twist DAs and are induced by the QCD EOM which relate twist-4 DAs 
to each other and to twist-2 and -3 DAs. The mass corrections to vector meson DAs 
have been calculated to twist-4 accuracy in Ref. [I6j. Quark-mass corrections to the 
evolution of DAs under a change of the renormalisation scale have so far only been 
calculated for twist-3 DAs [T7]. 

We shall study all these effects in this paper. The corresponding analysis of light pseudoscalar 
meson DAs, vr and K, can be found in Refs. [9l [HI [22] . 

In addition, in this work we present a new analysis of the parameters of the twist-4 DAs 
of the p meson. This update is long overdue: the "standard" values for these parameters 
can be traced back to a nearly 20-years-old work, Ref. [23] (see also [16]), and are in fact 
crude estimates obtained by dividing the leading QCD contribution to the relevant correlation 
functions by the typical hadronic scale. In this work we present, for the first time, a complete 
treatment of the twist-4 matrix elements within the QCD sum rule method. In particular, 
we resolve a puzzling discrepancy between the estimate of a certain next-to-leading order 
(NLO) parameter, in conformal spin, in Ref. [16] and in the renormalon model |13j . 

The presentation is organised as follows: in Sec. [2] we introduce notations and shortly 
review twist-2 and -3 DAs. In Sec. [3] we introduce the complete set of chiral-even and in 
Sec. [H chiral-odd twist-4 DAs. The conformal expansions of all DAs are worked out to NLO 
accuracy in conformal spin and reduced to a minimum number of non-perturbative hadronic 
parameters by solving the EOM constraints. In Sec. [5l we present models for these DAs, 
based on the calculation of the hadronic parameters from QCD sum rules. We summarise 
and conclude in Sec. [6] Details of the QCD sum rule calculations are given in the appendices. 

2 General Framework 
2.1 Kinematics and Notations 

Light-cone meson DAs are defined in terms of matrix elements of non-local light-ray operators 
extended along a certain light-like direction 2;^, 2;^ = 0, and sandwiched between the vacuum 
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and the meson state. Following Ref. [H], we adopt the generic notations 



(2.1) 



and 



M 



a] 



(2.2) 



for two- and three-particle DAs, respectively. The superscript A denotes the polarisation of 
the vector meson: A =|| (_L) for longitudinal (transverse) polarisation. The first subscript t = 
2, 3, 4 stands for the twist; the second, M = p, K*, . . ., specifies the meson. For definiteness, 
we will write most expressions for K* mesons, i.e. sq bound states with q = u,d. Whenever 
relevant, we will include quark- mass corrections in the form ± rriq, which allows one 



to obtain the results for mesons by rriq 



The variable u in the definition of two- 



particle DAs always refers to the momentum fraction carried by the quark, u = Ug, whereas 
u = 1 — u = Uq is the antiquark momentum fraction. The set of variables in the three-particle 
DAs, a = {ai, 0^2, 03} = {cts, C(q, ctg}, corresponds to the momentum fractions carried by the 
quark, antiquark and gluon, respectively. 

To facilitate the light-cone expansion, it is convenient to use light-like vectors and 
instead of the meson's four-momentum and the coordinate x^: 



m 



K* 



xP 



x^m 



K* 



X, 



X m 



2 ^ 

K* 



1 

-p,- + Oix' 
2 zp 



m 



K* 



pz 



(2.3) 



The meson's polarization vector e*-^-* can be decomposed into projections onto the two light- 
like vectors and the orthogonal plane as follows: 



pz 



e^^^p (X) e^'^^z j m 



pz 



pz 



2pz 



(2.4) 



We also need the projector g^^^ onto the directions orthogonal to p and z, 



9f,u {P,,Z^ + PuZf,) , 

pz 



(2.5) 



and will often use the notations 



a, = a^z^, bp = hf,p^ (2.6) 

for arbitrary four- vectors and h^. 

The dual gluon field strength tensor is defined as G^y = ^e^^p^G^"^. Our convention for 
the covariant derivative is = ~ ig^fi- Sometimes, a different convention for the sign 
of g is used in the literature: = + W^^. The sign of g is relevant for the parameters 
of three-particle DAs. 
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2.2 Conformal Expansion 

A convenient tool to study DAs is provided by conformal expansion, see Ref. [H] for a 
review. The underlying idea is similar to partial-wave decomposition in quantum mechanics 
and allows one to separate transverse and longitudinal variables in the Bethe-Salpeter wave- 
function. The dependence on transverse coordinates is formulated as scale dependence of 
the relevant operators and is governed by renormalisation-group equations, the dependence 
on the longitudinal momentum fractions is described in terms of irreducible representations 
of the corresponding symmetry group, the collinear conformal group SL(2,]R). The main 
rationale behind using the conformal expansion in the present context is that the EOM 
always relate contributions of the same spin; a truncation of the conformal expansion to a 
certain order is, therefore, consistent with the EOM. 

To construct the conformal expansion for an arbitrary multi-particle distribution, one first 
has to decompose each constituent field into components with fixed Lorentz-spin projection 
onto the light-cone. Each such component has conformal spin 

J = lil + s), 

where / is the canonical dimension and s the (Lorentz-) spin projection. In particular, 
1 = 3/2 for quarks and / = 2 for gluons. A quark field is decomposed as ip^ = A+ip 
and ip- = A-i' with spin projection operators A+ = 'jp'jz/ {^pz) and A_ = 'jz'Jp/ {2pz), 
corresponding to s = +1/2 and s = —1/2, respectively. For the gluon field strength there 
are three possibilities: Gz± corresponds to s = +1, Gp± to s = — 1, and both G±± and Gzp 
correspond to s = 0. Multi-particle states built of fields with definite Lorentz-spin projection 
can be expanded in irreducible representations of SL(2,M) with increasing conformal spin. 
The explicit expression for the DA of an m-particle state with the lowest possible conformal 
spin j = ji + . . . + jm, the so-called asymptotic DA, is given by [9J 

(pas{ai, aa, ■ ■ ■ , Om) = ^^i]^ ■ ■ ■ T{2j") ^i"'^^^"'^ ■ ■ ■ "m""^- (2-7) 

Multi-particle irreducible representations with higher spin j + n, n = 1, 2, . . ., are given by 
polynomials of m variables (with the constraint Ofc = 1 ), which are orthogonal over 

the weight function (12. 7p . For the two-particle DAs these are Gegenbauer polynomials, a 
convenient basis of orthogonal conformal polynomials for the three-particle DAs is given in 
App. A of Ref. [m (see also [21]). 

The anomalous dimensions of higher conformal amplitudes do, generally, increase loga- 
rithmically with the conformal spin, but a complete analysis of twist-4 anomalous dimensions 
is still lacking. It follows that all DAs approach their asymptotic form in the asymptotic 
limit as — > 0, i.e. at the scale /i oo. For practical applications one usually assumes that 
the conformal expansion is converging fast enough, so that a truncation after the few first 
terms is sufficient. The renormalon model of Refs. fi2[ [T3] presents an attempt to test this 
assumption, by giving an upper bound for possible higher-spin contributions. 
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2.3 Twist-2 Distributions 



The twist-2 DAs ^l.^, of K* mesons are defined in terms of the following matrix elements 
of non-local operators = 2u — 1) [15j: 



{0\q{xM-x)\K*{P,X)) = fl^rriK* l^P„, / due'^''^ 



e^^'>x 
Px 



{2.i 



\q{x)a,,s{-x)\K\P,\)) = 

^fk^[iel'^P.-ei'^P,)j\ue^^^^ 



{Px 



2 / ciMe^«^^ 



' II 1 1 



^ due'^''^ - + . . .| . (2.9) 



The above relations also include twist-3 and -4 two-particle DAs. The dots stand for further 
terms in which are of twist 5 or higher. The normalisation of all these DAs is given by 



(2.10) 



du (j){u) = 1 . 

The conformal expansion goes in terms of Gegenbauer polynomials: 

,fi) = 6uu\l + f2 ai'^(/i)Cf 2(2w -1)1. 



(2.11) 



71=1 



In this paper, we include terms up to NLO in conformal spin, i.e. truncate after n = 2. The 
dependence of the Gegenbauer moments a„ on the renormalisation-scale fi has been reviewed 
in Ref. [T7|, together with the numerical values of a„ and the decay constants /y"*"; these 
values are given in Sec. O 



2.4 Twist-3 Distributions 

To twist-3 accuracy, there is a total of four two-particle DAs and three three-particle DAs. 
Two of the former, (p^.^* and ^f./^*, have already been defined in the previous subsection. 
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The other two are given by 



2 

"1 



(2.12) 



mz)s{-z)\K*{P,X)) = -z/^,(eWz)m^, / due'^P'4^,{u) (2.13) 

Jo 

with the normahsation 

rd.c!(")=i-C^^. (2.14) 

Jo fp>, ^K* 

The three-particle DAs are given by: 

{Q\qiz)gGp,ivz)^,^,si~z)\K*iP,X)) = /imx.(pz)H'^$|,,.(t;,p^) + . . . , 

{Q\q{z)gGp,ivz)tj,s{-z)\K*iP,\)) = /imx*(p;^)H>S;i^.(^>P^) + • • • , 
{0\qiz)gG,pivz)a,psi-z)\K*iP,\)) = /^,m^,(e(^)^)(p^)$3^;,.(t;,p^) , (2.15) 
where the dots denote terms of higher twist and we use the short-hand notation 

J^{v,pz) = y"pae"*P'("2""i+^"^)j^(a) (2.16) 

with J^{a) a three-particle DA. a is the set of parton momentum fractions a = {ai, a^-, 03} 
and the integration measure Da is defined as 

y T)a = J daida2da3 6 (^1 — ''^^ . (2-17) 

The above DAs are not independent of each other, and their mutual interrelations have 
been unravelled in Ref. |15j, including quark- mass corrections. Explicit expressions for the 
conformal expansion are given in Ref. [17], together with the /i-dependence of hadronic 
parameters. Numerical values are given in Sec. [H 



3 Chiral-Even Twist-4 Distributions 

In this section we derive expressions for the chiral-even two- and three-particle twist-4 DAs 
of the K* to NLO in the conformal expansion. The corresponding expressions for the p 
were obtained in Ref. [16]. In this paper we include also G-parity- violating and explicit 
quark-mass corrections. We compare the resulting DAs obtained in conformal expansion 
with those following from the renormalon model developed in Ref. [13] . Numerical estimates 
of the hadronic input parameters and the resulting DAs are discussed in Sec. El 
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We start with the three-particle distributions. The analysis closely follows that of Ref . [16] . 
There are four chiral-even K* three-particle DAs of twist 4, defined as [16] 

{0\qiz)^,^,gG^pivz)si-z)\K*iP,\)) = 

, ^^^^^ 



pz 



{0\q{z)zj^gG^p{vz)s{-z)\K*{P,X)) = 



(3.1) 



/ \ ^^^^ z 

Pt^ (e!£p/3 - e^/3Pa) fLmK*'l>l.j^,{v,pz) + {pag^^ - Ppga/,) " — fi*m%.<^lj^.{v , pz) 



+ Pm {PaZfi - PisZa) fK*rnj^,^l.j^,{v,pz) + 



(3.2) 



the dots denote terms of twist 5 and higher. 

\&4.j^* and correspond to the light-cone projection "jzCzp which picks up the s = 

+1/2 components of both quark and antiquark field and the s = component of the gluon 
field. The conformal expansion then reads: 



^4;i^.(^) = 120aia2a3[V'! + ^i(ai - "2) + V'I(3«3 -!) + •••], 
5lK*(a) = 120aia2a3[V'o + ^i(«i - "2) + ^|(3a3 - 1) + . . .]. 



(3.3) 



G-parity implies that, for the p and meson, ipQ='il'2—'^i— 0; whereas for the K* meson 
ipl, ipl and ipf are 0{'ms — m^). 

In turn, the DAs $4;x* and ^^-k* correspond to the light-cone projection 'j±Gz_i, which 
is a mixture of different quark-spin states with Sg = +1/2, Sg = —1/2 and Sq = —l/2,Sq = 
+1/2, respectively. In both cases s = +1 for the gluon. We separate the different quark-spin 
projections by introducing the auxiliary amplitudes and , defined as 

{0\q{z)gG,,{vz)jzlal,lps{-z)\K*{P,X)) = flm%.{e^'h){p,g^, - p,g^W^{v,pz) , 



{0\q{z)gG,,{vz)jpj^j,jX-z)\K*{P,X)) = f^,m^K*ie^^^ z){p,gt ' Pug^J^" iv,pz) . 



The distributions $4.^^* and $4.^^* are then given by 

1 



a) 



a 



(3.4) 



(3.5) 



1 In the notation of Ref. [16], $|{.^. = $, *|.^^. = = 4>, ^[.j^, = ^ 
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$1^^ and have a regular expansion in terms of conformal polynomials: 



= 60a2«3 



+ 0i' (as - 3a2) + 02^ ( "3 - T^tti 



For the p and meson, G-parity implies 



at 



so 



that (f)]^ = (f)l^^ For K*, we write 

0P = 0| + ^1 



at 



(3.6) 



(3.7) 



(3.^ 



where the are the G-parity- violating corrections. Using (13.51) . we readily derive the fol- 
lowing expressions: 



^lK*i^) = SOttg <^ 0{i(l - as) + (pl [^3(1 - as) - Gaiaa] 



+ 



"3(1 - "3) - 2(^1 + "2) 



(ai - a2 



^0 + asOl + - (Sag - 3)^^ 



= SOag |6'o(l - as) + [asi^l - as) - 6aia2] 



as(l - as) - -{al + al) 



- (ai - a2) 



t>o + a30i + - (5as - 3)02 



In addition, we introduce one more three-particle DA (a) [13] : 

(0 |g(2;)7c,^L'^G^i.(f2;)s(-z)| K*{P, A)) = fl^m^K^PaPu El.j^,{v,pz) + 

pz 



(3.9) 



(3.10) 



The Lorentz structure PaPu is the only one relevant at twist 4. Because of the EOM D^G^ 



9 So 'y/sq, where the summation goes over all light flavors, S^j.^, can be viewed 



as 



describing either a quark- ant iquark-gluon or a specific four-quark Fock-state of the K*, with 
the quark-antiquark pair in a colour-octet state and at the same space-time point. The 
conformal expansion of starts with the conformal spin j = 4 and reads 



S||.^,(a) = 840aia2a3 



(3.11) 



where Q is dimensionless. The dots stand for terms with higher conformal spin j = 5, 6, . . ., 
which are beyond our accuracy, sf was not considered in Ref. [16] because is G-odd and 
the first G-even term only occurs at the next order in the conformal expansion, for j = 5. 



^This implies, in particular, that only one of the DAs ^4-p(^) and $4;p(0) is dynamically independent. 
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Eqs. (13. 3p . (13. 9p and (13.111) represent the most general parametrization of the chiral- 
even twist-4 DAs to NLO in the conformal-spin expansion and involve 13 non-perturbative 
parameters. Not all of them are independent, though. In the following, we shall establish 
their mutual relations and also express all leading order (LO) and the G-even NLO expansion 
coefficients in terms of matrix elements of local operators. 

Except for the asymptotic three-particle DAs correspond to contributions of the 

lowest conformal spin j = js + jq + jg = 3. The parameters ipl, ipoy 0o multiplying 
the asymptotic DAs can be expressed in terms of local matrix elements as 



{0\qgG^pz^,s\K*{P,X)) = /]! m^.^^. je^ (^P^P, -^ml,g^,^ - (a ^ P) 



+ I fim\*4K* {e'^a^ap,. - ef^'ga^^ ■ (3.13) 

Here we adopt the generic notation that ^ are G-conserving and k G-breaking parameters. 
(^3, K3 are twist-3 and (4, K4 twist-4 parameters. 

Taking the local limit of Eqs. (13.11) and (13.21) . and comparing with the above definitions, 
one obtains 

00 = CIk* + 2 (Ik* ' ^0 = ~2 '^Ik* + 2 '^Ik* ' 

ii 2 ii 1 ii ~|| 2 ii 1 ii 
Vo= 3 4k* + 3 4k* , Vo= ^ Qk* + 3 Qk* ■ (3-14) 

The results for ^jj and ipQ agree with those obtained in Ref. [16], the others are new. Note 
that the "twist-4" DAs receive contributions from both twist-3 and -4 operators. This is due 
to the fact that the standard counting of twist in terms of "good" and "bad" components 
introduced in Ref. [25] differs from the definition of twist as "dimension minus spin" of an 
operator. See also the discussion in Sec. 2.2 of Ref. [15] . 

What about the scale-dependence of these parameters? The relevant local twist-4 op- 
erator mixes with operators of lower twist for 7^ 0. Neglecting 0{m1) corrections, the 
mixing is given by 

{qial,gG,.sf = {qi^l,gG,^sf^ ("i _ ^ ^ In 4) + ^ - 1^ 4 ■ (3-15) 

The matrix element of the derivative operator on the right-hand side vanishes for vector 
mesons, so that Cjx* renormalises multiplicatively even for rris 7^ 0. Resumming the loga- 
rithm, to LO accuracy, one has 

CUf^') = L''^^''"k!K*{l^l) , (3.16) 
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with L = as{iJ^)/as{ixl). 

The scale dependence of most easily be derived by observing that this parameter 

is related to af and quark masses by the QCD EOM [19] : 



II 3 II f^* ms - ma ml - ml , ^ 

J = J_:q^^ ^+ ' ^ 3.17 

Taking into account the known scale dependence of a|, and m^^g, one obtains 
AkAi^') = AM) - (L3^/(^^«)-l)a!(/.g) 

_ /^16/(3A,) _ /i*(/^o) K -mg](/ig) ^ _ ^\ K -"-^^ll^o) 



.(3.18) 



To NLO in conformal spin, the discussion becomes more involved. As explained in Ref . [H] , 
for massless quarks the corresponding contributions can be expressed in terms of matrix 
elements of the three possible G-parity-even local quark-antiquark-gluon operators of twist 
4. These three operators are not independent, however, but related by the QCD EOM. One 
is left with only one new non-perturbative parameter, Zj\j^» jf| which can be defined as 

- r ~ 1 4 

(0|g Y^^.gG^f^^ 7575S - -{id^qgG^frifr^^s\K\P, A)) + (/i ^ z/) = 

= 2fl.m\, [efP, + eWp^) . (3.19) 
The scale dependence of , for massless quarks, is given by 



i5|^,(/i2)=L^°/'3»i3l^.(/i, 



For massive quarks, the twist-4 operator mixes with operators of lower twist. These lower- 
twist operators have the same Dirac structure as in (13.151) . but additional derivatives acting 
on the fields. This means that, in terms of DAs, mixes with V'ii^*) Eq. (I2.13p . which in 

turn mixes with the three-particle twist-3 DA $3.^, , Eq. (I2.15p . which itself mixes with twist- 
2 DAs fTT]. As the numerical impact of the admixture of m^ times lower-twist parameters is 
negligible for all cases investigated so far (twist-3 and -4 pseudoscalar parameters [11] and 
twist-3 vector parameters [13 )? we refrain from working out these relations. 

For massive quarks, one has to distinguish between G-parity-conserving and G-parity- 
breaking contributions. G-parity-conserving corrections do not involve new operators, and 
the difference to the massless case is mainly due to corrections proportional to the meson 
mass. This case is described in detail in Refs. [T31[TB]. Here we just quote the results obtained 
in Ref. [16]: 

,11 _ 1 II 5 II 3 ^11 1 II 7^11 



^In the notation of Ref. [12], ^\k- — ^'t'^t- 
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3 
4 



1 

- UJ. 



3K' 



+ 7uj 



21 
T 



7"- 



(3.20) 



where the terms in a|, C|x*' '-^3i<'* ^^"^ ^lic* mass corrections, wl^. and oolx* are twist-3 
parameters and have been discussed in Ref. [17]. 

The G-parity-breaking contributions, on the other hand, involve a different set of local 
operators and in particular 

^=u,d,s 

which determines the normalization and the leading conformal spin contribution to the DA 
defined in Eq. fl3.10p . Hence, a complete treatment of G-parity-breaking corrections 
to twist-4 DAs requires also the inclusion of It is beyond the scope of this paper 

to work out the corresponding relations between the matrix elements of local operators 
and expansion coefficients. Instead, we adopt a different approach and estimate G-parity- 
breaking corrections of conformal spin j = 4 using the renormalon model of Ref. [13] . The 
general idea of this technique is to estimate matrix elements of "genuine" twist-4 operators by 
the quadratically divergent contributions that appear when the matrix elements are defined 
using a hard UV cut-off, see Ref. [13] for details and further references. In this way, three- 
particle twist-4 DAs can be expressed in terms of the leading- twist DA 



,R 
4;K* 



(«) = - 



,R 
4;K* 




"2</>2:X*("l) ai02-ft:*(«2) 



1 — «! 



1-02 



(3.21) 



where, in contrast to Ref. [13], we do not assume that 0|.^.(m) is symmetric under the 
exchange u ^ u. 
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The expressions in fl3.2ip do not rely on conformal expansion and contain the contribu- 
tions of all conformal partial waves. Projecting onto the contributions with the lowest spin 
j = 3, 4, we obtain 



^0 = 


0, 


^1 = 






7 li.ll 




-c" 


7113 


7 11.11 


71!. R 




00 = 




0^ = 










0, 




7 11.11 


^|'^ = 


7 11.11 


so — 


1 11.11 











(3.22) 



These results have to be compared with those in Eqs. ( 13.14p . ( ]3.20p in the limit that all con- 
tributions from twist-2 and twist-3 parameters are set to 0. It follows that in the renormalon 
model 

_ll 1 

uJai^, = -- , (3.23) 



1 

^4K' - -g , 

which has the same sign as, but is larger than the result from a QCD sum rule calculation 



see Sec. Also note that in the renormalon model = = 0, in contrast to fl3.14p . This 



is due to the fact that the contribution in in (13.141) is obtained as the matrix element 
of the operator (13.131) which vanishes by the EOM (up to a total derivative), see (13.171) . 
Therefore, against appearances, this contribution has to be interpreted as "kinematic" mass 
correction induced by the non- vanishing i^'*-meson mass rather than a genuine twist-4 effect. 
The complete expressions for the G-odd parameters V'l, V"!; ^1,2 ^0 ^^^^ contain mass 
corrections in terms of lower- twist parameters, whose determination is, as said before, beyond 
the scope of this paper. 

We are now in a position to derive expressions for the chiral-even two-particle DAs of 
twist 4, and ipl.j^*, which are defined in Eq. (12.81) . From the operator relations collected 
in App. [X]we obtain: 



u 



du 



dai 



da2 — 

"3 



2$L,(«) + ^ 



4::K* 



.a 



rUg f^, d 



rriK' 

2^ ( 4;K* 
d 



K* 



(3.24) 



u) - 



dv 



4— / dai 

du 



da2 — («! - ^2 - 



On 



2$" («) + ^i 



4;_ft:* 



a 



12 



rriK* /|.. du 



+ ^^^^^3^2^" . (3.25) 



The latter equation has to be integrated with the boundary conditions 04-;^. (0) = = 
0|.^,(1) which imphes the relation 03.171) between a| and The boundary condition 

arises from the conversion of the matrix element of Eq. ( ]A.1|) for the K* into a relation for 
(t>\-K*- This derivation of fl3.17p is equivalent to that given in Ref. [T9j . 

iA-K* corresponds to the projection s = — | for both quark and antiquark and hence, in 

1 /2 

the absence of quark-mass corrections in ± m^, has an expansion in terms of C„'"(0. It 
can be split into contributions from genuine twist-4 parameters which are defined in terms 
of local twist-4 operators and kinematic Wandzura-Wilczek-type and mass corrections: 

(«) = ^1?* («) + V'I;x*^(m) • (3.26) 

V'l'^.^ contains corrections explicitly proportional to rru ± m^, of which we only keep the 
leading terms in (m^ ± m^)^, but neglect higher powerso We then find 

^28 4 ~ T ~ 16 / 



+ 6^^^^#^|e + a^^(3e^-l) + a^^e(5e^-3) + ^.3V(3e^-l) 



+^ ^3Ve(5e' - 3) - ^ (35^^ - 30^2 + 3) ^ . (3.27) 



04.^*, on the other hand, has no regular conformal expansion and contains logarithms even 
in the chiral limit. We solve the integral relation (13.251) by substituting all DAs on the 
right-hand side by their conformal expansion to NLO, implementing the boundary condition 
04.i<-.(O) = = 04.^.(1) by eliminating in favour of af according to (13.171) . and dropping 
terms in (m^ ± m^)" with n > 1. Like with V'l ii-*) ^6 distinguish between genuine twist-4 
and mass corrections and write 



M=0Kt(«) + 0S'^M- (3.28) 



The terms in (m^ ± TOq)^ actually diverge for u — > 0, 1, which is however largely irrelevant for phe- 
nomenological applications. 
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We then find: 



84(^V i I (21 - 13C^) + M^(10 - 15m + 6u^) lnu + #(10 - 15u + Gm^) Inw 
8 

+ 80^'^ \ u^{2 -u)\nu- #(2 -u)\nu-l (3^^ - H) i , 



"^10 (7 9 6 ~ 4 '^^^*) ^2^^(0| 

+ 2 I -2ai - ^ CL. + 3u;S^. | 1 1 «f2 (21 - 13^') 

W{10 - 15u + 6w^) lnu + u^{10- 15u + 6fZ^) In u} 
+ 4 |al - y 4k* I |^^(2 -u)\nu- #(2 - fZ) InfZ - ^ (3^^ - 11) | 

+ UlllJIh. IfL Quu (2(3 + 16ai) + ^(4^, - ai)C!^\0 



rriK* /II. L 



+ 24 1h±J^ 1(1 _ 3ajL + 60^)1^2 In + (1 + 3ajL + 60^)1^2 In u] 

+ ^^^^ 4^ 6. J- f 10.3V + f ai] cni) 



Ik 

+ 20 + i a- - 1 (0 + A3V + ^ C^f (0 

■'- -L ^ , \ /-v3/2/>.x , 2 _L 3/2, 



315 - 21 J (0 + ^ A3vcr(e) 

+ Hh^:^ if- {{5u^ - 23 - 54a^ - 108a^) Inu 

rriK' /II, 

- {5u^ - 23 + 54a]^ - 108a^) hiu} . (3.29) 
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Recall that u = 1 — u and ^ = 2u — 1. The DAs for K* = (qs) mesons are obtained by 
replacing uhj 1—u. Both (13.261) and fl3.28p agree, for the p meson, with the results obtained 
in Ref. [16]. 

The above expressions provide a self-consistent model of the twist-4 DAs which includes 
the first three terms of the conformal expansion. As mentioned above, one shortcoming of 
the model is that G-parity-breaking terms in tpl and Of 2 are not known exactly, but only in 
the renormalon model, Eq. fl3.22p . which misses meson mass corrections. Numerically, the 
neglected parameters may be of the same size as the included ones. 

An estimate of the size of higher orders in the conformal expansion can be obtained using 



T4 



the full renormalon model for ^4.;^ 

2 



;II:T4,R/ X 



^ll,T4,R 



c 



iK 



and In this case, one has 



— u 



u) 



~>AK* 





u + u + (m — w) In 



u — V 



+ 



dv 



u + + {v — u)\ia 



u 



(3.30) 



As explained in Ref. [13], the renormalon model does not take into account the damping of 
higher conformal-spin contributions by the increasing anomalous dimensions and, therefore, 
only provides an upper bound for their contribution. The most important effect of these 
corrections is to significantly enhance the end-point behaviour of higher-twist DAs in some 
cases, which can be important in phenomeno logical applications, for instance Z) — > (tt, K) 
form factors [26] . 



4 Chiral-Odd Distributions 

The analysis of chiral-odd DAs proceeds along similar lines and, except for the inclusion 
of G-odd contributions, replicates that performed in Ref. [1^. We first define the relevant 
three-particle DAs: 

{mz)a^p9G,u{vz)s{~z)\K*{P,\)) 

= fk*i^h [PaPi,9p^ - ppPi^gL - Po^PuOp^, + ppPugil,] ^i-K* {v, pz) 

, 2 r (^) -L W -L (A) ± , (A) ±i^±(2)/ N 

fK*mK* I (A) (A) (A) , (A) i^±(3)/ x 

+ [PaPf,e\fiZ^ - PpPi,e\;^z^ - PaPveYpZf, + ppp^ey^Zf,\^^.y.{v,pz) 



15 



/^*m^. (A) (A) (A) , (A) i^-L(4)/ X , |.^^ 

+ [PaP^,e\lzp - ppp^,e\lza - PaPuey^zp + ppPuey^Za\^^^j^',{v,pz) + ... , (4.1) 

pz 

m{z)9G,u{vz)s{-z)\K*{P,\)) = tf^,m'j,4e^^y-e^^lp,]^i^^,{v,pz) + ... , 
{0\q{z)tgG,,ivzh,s{-z)\K*{P, A)) = z/^.m^^e^k - e^llp,]^J,,{v,P^) + ■■■, 



{0\q{z)ai,^gDaGaf^{vz)s{-z)\K*{P,X)) = if^.mj^, e^±lPiy - e''_^lPi^ Pi3^t-K*{^^Pz) + 



(4.2) 



The twist-3 DA already mentioned in Sec. [21 the twist-4 DAs are related to those 

defined in Ref. [TB] by = t/^^ ^ix* — "^i-K* = S^^. was first introduced in 



Ref. [13]. As usual, the dots denote terms of higher twist. 

As the matrix element in (14. ip is G-odd, the DAs ^^^5 limit of equal mass 

quarks, antisymmetric under the exchange of ai and a2, whereas "^^-k* ^ind ^1'^^. are sym- 
metric. In order to resolve the conformal structure of useful to exploit the fact 
that (Jfj^ulb is not independent of a^i,, and to define the "dual" matrix element 

m{z)i(yapl^gG^,{vz)s{-z)\K*) = r.h.s. of ^ with . (4.3) 

One easily finds 

^± _ ^± ^±(1) _ _^±(3) S±(2) _ _^±(4) 

$^(3) _ g±(4) _ _^±(2) . . .^ 

and $^^1 correspond to the Lorentz spin projection s = +1/2 for both quark fields 
and s = for the gluon. Hence the conformal expansion reads, to NLO: 

^iiJ-U") = 12Oaia2a3[0o + - "2) + 0^(3^3 - 1)] , 

(a) = ^t^l{a) = 12Oaia2«3[0o + - "2) + 0^(3^3 - 1)] • (4.5) 

Here ^^re G- violating and 0]*- is G-conserving; the same applies to 0^. 

The conformal expansion of S^^^, starts at j = 4 and reads, to NLO accuracy: 

^i,K' (a) = 8A0aia2al^^ . (4.6) 

The remaining DAs do not correspond to fixed values of the Lorentz-spin projection s. 
Like in the chiral-even case, we introduce auxiliary amplitudes with a regular conformal 
expansion: 

{0\qizh,^pgG,,ivz)s{-z)\K*iP,\)) = tf^,mj,.{pz)[e^iy-e^llp,]^^^{v,pz), 
{0\qizh,^,t^,gG,,ivz)si-z)\K*iP, A)) = tfk,mUp^)[e^iy - e^llp,]^^^{v,pz), (4.7) 
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and, similarly, two more distributions "^^^ and by replacing 7^7^ ^p^z- The DAs in 
(14. ip and (14.21) are then given by: 



and have a regular expansion in terms of conformal polynomials, to wit: 



tPl^ + iPl^ias - 3ai) + ^\)\' ( as - 



(4.8) 



V^o^^ + 7/'jT(a3 - 3«2) + ^2' ( "3 - 



(4.9) 



Again G-parity ensures that for the p and mesons we write 

^^^=^i^Qi. ^f = ^.^-^^ (4.10) 

where the Q^- correspond to SU(3)-breaking corrections that also violate G-parity. Introduc- 
ing an analogous decomposition of and in terms of ip^- and ^ we then find 

'^^-kA^L) = SOftg {i^oil - as) + [03(1 - as) - Gaiaa] 

3 



+ ^2 



"3(1 - "3) - 2("i + "2) 



(ai - 02) 



+ ag^i + 2 (5«3 - 3)^2 



'^i-K'i^ = SOog |^c)"(l - "3) + i^i ["3(1 - "3) - 6aia2] 



+ ^2 



03(1 -«3) - -(01+^2) 



(ai - 02) 



+ 03^1^ + 2 (5«3 - 3)^2 



= -30a^ {^0^(1 - ^3) + ei [«3(1 - «3) - 6aia2] 



+ ^0 



03(1 - ^3) - -(«! + 



(ttl - ^2) 



w + oi^i^i + 2 (^"3 - 3)V^^ 



$^(Jl(a) = 30a32 {^0^(1 - a^) + 0^ [a3(l - as) - 6«ia2] 



+ ^0 



a3(l - ^3) - \{al + al) 



(ai - ^2) 



(4.11) 



Our next task is to relate ifjf- and 9^ to local matrix elements. To LO in the conformal 
expansion, the relevant matrix elements are [16] 



mgG^,^s\K\P, A)) = ifkml.CiK* U'^Pp - Pc 
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{0\q9G^,afs,s\K*iP, A)) = /^m^. jl^^. (e^P^ + ef^P^) + k^^. (e^P^ - ef^P^) | , 

(4.12) 

where again C, denotes G-parity-conserving parameters and k G-parity-breaking ones. The 
twist-3 parameter Kg^, was investigated in Ref. [l7j. Taking the local limit of (14. ip and (14. 2 p 
and comparing with the above, one finds 



^0"- = 








~ Cik* 5 




00^ = 






K 




1 

K 

3 ' 


0^ = 




- V 


0^ 




1 



(4.13) 



Whereas Cik*; Cik* ^'^d k^j^, are independent parameters, k^^, depends on quark masses 
and aj^ by a relation that is analogous to Eq. (13.171) and was obtained in Ref. [21] : 



- 10 + /i, 12m^. Ami. ' ^ ^ ^ 

Like with the scale-dependence of Kf^-. follows from that of the parameters on the 

right-hand side and is given by 



t4A-.(f') t„.,„„ , ^_ ,^„, , "'/i.(Aii) 12mK. 

- (L«/«. - 1) 1'"--'"'"''°' (4.15) 

with L = as(/i^)/a;s(/io)- In the limit of zero quark mass, the parameters Cj, Cl renormalise 
multiplicatively ^23j : 

(cr + ci) (/i^) = (cr + cr) (z^^), 7+ = 3^^^ - ^ c^, 

(cr - cr) (/i^) = (cj - cr) (i^d, 7- = - ^c,. (4.16) 

Again, these simple scaling relations will receive corrections from terms proportional to the 
quark masses. These corrections are unknown, but based on the experience with pseudoscalar 
twist-4 matrix elements [13] we do not expect them to be large. 

The calculation of the NLO G-even corrections is pretty involved and presented in detail 
in App. B of Ref. [IB]. The upshot is that the six coefficients ipi2y '^12^ ^1 ^1 involve 
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three additional genuine twist-4 parameters {{Q^^'^'^^)): 



^2^ 



-^4 + ^.3V-^((Q«)) + ^((Q^^^)) + ^((Q^^^)) 

44-^^3V + |((Q«)) + ^((Q^^^)) + ^((Q^^^)) 



^2 = «2 - ^, 



22 



12 



'^3ii'* 



55 



11 



(4.17) 



The precise definition of ((Q*-*^)) is given in Ref. [I6j. ijJ^x* is a twist-3 parameter and defined 
in Ref. [17]. Existing numerical determinations of these parameters from QCD sum rules 
are far from being precise, so we decide to estimate them using the renormalon model of 
Ref. [13] instead. The model entails the following expressions of the three-particle twist-4 
DAs in terms of cj)2-K* '■ 



^±(l),R/ N ^±(3),R/ X 



«20^K.(ai) ai0^x*(a2) 



;i -"2) 



^±(2),R/ X 



a) 



.V 



1 -ai) 



(1 - "2) 



1 — «! 



+ 



1-^2 



a2 



1 — ai 



1-^2 



(4.18) 



which is similar to the results for the chiral-even DAs, Eq. (13.211) . The model implies the 
following estimates for the G-conserving twist-4 parameters in (14.13^ and ( 14. 17^ : 



c 



,R 



< 



AK* ) 



{i)\\R 



10 



-C 



4:K* 5 



{{Q^'^r = 0. (4.19) 
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The G-parity-violating parameters are given by: 



^0^'^ = , 



_,R 



^1^ 



±,R 



g S4X*"l 5 SO 



21 M ± 2^ 
20 ' '^'^ 



,R 



,R 



20 ^^^^^ ' 



(4.20) 

As with chiral-even distributions, the above results provide an estimate only for the genuine 
twist-4 contributions, but miss any mass corrections in terms of lower-twist parameters. This 
is also the reason why (pQ, 6q, 0j} and ^^j" vanish in the renormalon model, in contrast to 
Eq. 

Let us now turn to the two-particle twist-4 DAs 0^^. and V'i^^* defined in Eq. (12.91) . 
From the operator relations given in App. |X]we obtain: 

d n f 1 

-2 — / dai / (ia2 < ^ («! - ^2 - 0^ix*(^) 



du 



"3 ^ 



+ 2 



03V(n), (4.21) 



d r , r , 4 f a,-a,-^ 

;— / dai / da2 — < 

Jo "3 I "3 



du 



'0 

+ ^03.;^*(m). 



_(3) 



(4.22) 



The boundary condition for (p^.j^* is (pi-^'i^) = = 0^^.(1), which implies the relation 
(14.141) . In the renormalon model, one obtains exact expressions for these DAs [T3] : 



..T4,R, 



U) = 8(4^. 



dv [ u + {u — v)\n 



u 



/ dv i u + (v — u)\n ) 

Ju \ V J 



(4.23) 



Like the chiral-even DA ipl.j^*, ipi-x* corresponds to the projection s = — | for both 
quark and antiquark and hence, in the absence of quark-mass corrections in ± rriq, has an 
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expansion in terms of Cn'^{^). The full i^^.x* contains corrections explicitly proportional to 
iTLs ± iriq, of which we only keep the leading term in (m^ ± . To NLO in the conformal 
expansion, (14.211) yields: 

^iMu) = 1 + (^124^. - ^ ai^ cl'\i) + [-1 + ^-4- 10 {ci,. + Hk*]) cl'\o 



3A'' 



124^. --A^^, + -a^ + 5 



1^ {3 (1 + 6a|) + ?>a\cl/\0 + 5 (4CL* - 3a|) ^^'/'(O 



6 



+ 6 



MM ( 9a" + 10^02 



I ^1 — 3af + 6a|j m Inw — ^1 + 3af + 6a|j m Inuj . (4.24) 



Recall that u = 1 — u and C, = 2u — 1. The above expression refers to a K* = (sq) meson; 
for K* = (qs), one has to replace uhj 1 — u. 

The explicit formula for 0^^. from (I4.22p is very long and complicated, so we only give 



the result for ± rriq 0: 



^ 25 ^ 3 45 15 ^ 30 ^ 5 ^ 2 ) 1 v,; 



* 35 ^ 60 / ^ 1575 ^ 



+ (5/?^^. -a^-200^) <j -4^3(2 -M)lnM + 4^3(2 -M)lnu + -M<(3^2- 11) 



+ (2<,.-??a2--^((Q«))-^((g(3))) 



11 



55 



11 
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X < u^{6u'^ - 15m + 10) liiM + M^(6u2 - 15m + 10) Inn - - uU (13^^ - 2l) 



(4.25) 



Both fl4.24p and fl4.25l) agree, for the p meson, with the results obtained in Ref. [ISj. The 
numerics of the above DAs will be discussed in the next section. 

5 Models for Distribution Amplitudes 

In this section we compile the numerical estimates of all necessary parameters and present 
explicit models of the twist-4 two-particle DAs introduced in Sees. [3] and |H The important 
point is that these DAs are related to three-particle ones by exact QCD EOM and have to 
be used together: this guarantees the consistency of the approximation. Our model thus 
introduces a minimum number of non-perturbative parameters, which are defined as matrix 
elements of certain local operators between the vacuum and the meson state, and which we 
estimate using QCD sum rules. More sophisticated models can be constructed in a systematic 
way by adding contributions of higher conformal partial waves when estimates of the relevant 
non-perturbative matrix elements will become available. 

Our approach involves the implicit assumption that the conformal partial wave expansion 
is well convergent. This can be justified rigorously at large scales, since the anomalous 
dimensions of all involved operators increase logarithmically with the conformal spin j, but 
is non-trivial at relatively low scales of order fi ~ (1^2) GeV which we choose as reference 
scale. An upper bound for the contribution of higher partial waves can be obtained from the 
renormalon model. 

Since orthogonal polynomials of high orders are rapidly oscillating functions, a trun- 
cated expansion in conformal partial waves is, almost necessarily, oscillatory as well. Such 
a behaviour is clearly unphysical, but this does not constitute a real problem since physical 
observables are given by convolution integrals of DAs with smooth coefficient functions. A 
classical example for this feature is the 77*-meson form factor, which is governed by the 
quantity 



where the coefficients are exactly the "reduced matrix elements" in the conformal expan- 
sion. The oscillating terms are averaged over and strongly suppressed. Stated otherwise: 
models of DAs should generally be understood as distributions (in the mathematical sense). 

All relevant numerical input parameters for our model DAs are given in Tabs. [1] and [21 
at the scale /i = 1 GeV, which is appropriate for QCD sum-rule results, and at the scale 
/i = 2 GeV, using the scaling relations given in Sees. [3] and HI to facilitate the comparison 
with future lattice determinations of these quantities. 

The parameters related to twist-2 matrix elements have been determined using various 
methods; see the discussion in Ref. [17]. Matrix elements of twist-3 operators were also 
discussed in Ref. p7]. Twist-4 matrix elements for the p were estimated a long time ago 
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P 


K* 







fx = 1 GeV 


/i = 2 GeV 


/i = 1 GeV 


/i = 2 GeV 


fl = 1 GeV 


/i = 2 GeV 


fl [MeV] 


216(3) 


216(3) 


220(5) 


220(5) 


215(5) 


215(5) 


[MeV] 


165(9) 


145(4) 


185(9) 


163(8) 


186(9) 


164(8) 


a\ 








0.03(2) 


0.02(2) 








ai 








0.04(3) 


0.03(3) 








al 


0.15(7) 


0.10(5) 


0.11(9) 


0.08(6) 


0.18(8) 


0.13(6) 


ai 


0.14(6) 


0.11(5) 


0.10(8) 


0.08(6) 


0.14(7) 


0.11(5) 




0.030(10) 


0.020(9) 


0.023(8) 


0.015(6) 


0.024(8) 


0.017(6) 










0.035(15) 


0.017(8) 










-0.09(3) 


-0.04(2) 


-0.07(3) 


-0.03(2) 


-0.045(15) 


-0.022(8) 










0.000(1) 


-0.001(2) 










0.15(5) 


0.09(3) 


0.10(4) 


0.06(3) 


0.09(3) 


0.06(2) 


^3V 








-0.008(4) 


-0.004(2) 
















0.003(3) 


-0.001(2) 










0.55(25) 


0.37(19) 


0.3(1) 


0.2(1) 


0.20(8) 


0.15(7) 










-0.025(20) 


-0.015(10) 









Table 1: Decay constants and twist-2 and -3 hadronic parameters at the scale = 1 GeV and 
scaled up to /i = 2 GeV. The sign of the twist-3 parameters corresponds to the sign convention for 
the strong coupling defined by the covariant derivative = — igA^f^; they change sign if g is 
fixed by = dfj, + Numbers taken from Ref. [T7], see also Ref. [8] for a detailed discussion 

of the decay constants. 



from QCD sum rules [HI [23l ETj. In this paper, we perform a complete reanalysis of these 
parameters and also include G-parity-breaking effects relevant for the K* and SU(3) breaking 
for the meson. The corresponding sum rules and plots are given in the appendices. 

For the chiral-even parameter (1^ we find Cip = 0.07 ± 0.03, which agrees with our older 
result dp = 0.15 ± 0.10 [T6] within errors. The change is due to updated input parameters. 
Another parameter, a;|^, was estimated, in Ref. [16], from a correlation function of currents 
with different chirality, by dividing the leading contribution (a dimension-5 power correction) 
by the typical hadronic scale. The result u;|^ = 0.1 ± 0.1 is a crude estimate. In this paper 
we obtain = —0.03 ± 0.01, from a careful analysis of various sum rules. This result is 
smaller than the previous estimate and negative, in agreement with the prediction based 
on the renormalon model fl3.23p . Th absolute size is smaller than the renormalon-model 
prediction, which is not significant, however, as the intrinsic renormalisation scale at which 
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p 


K* 





H = 1 GeV 


fi = 2 GeV 


/i = 1 GeV 


= 2 GeV 


/i = 1 GeV 


12 = 2 GeV 


d 


0.07(3) 


0.06(2) 


0.02(2) 


0.02(2) 


0.00(2) 


0.00(2) 




-0.03(1) 


-0.02(1) 


-0.02(1) 


-0.01(1) 


-0.02(1) 


-0.01(1) 

V / 


Ci 


-0.03(5) 


-0.02(3) 


-0.01(3) 


-0.01(2) 


-0.01(3) 


-0.01(2) 


Ci 


-0.08(5) 


-0.05(3) 


-0.05(4) 


-0.04(2) 


-0.03(4) 


-0.02(2) 










-0.025(5) 


-0.020(4) 








nix* 








0.013(5) 


0.011(5) 









Table 2: Twist-4 parameters at the scale /U = 1 GeV and scaled up to 2 GeV. Sign convention for 
the strong coupling g as for twist-3 parameters in Tab. [H 



the model is valid is not fixed. 

Another important result is that we find Ci" + Ci" 0- This parameter is usually set 
to zero, based on the observation that the leading contribution to the correlation function 
of the corresponding quark-quark-gluon operator with the electromagnetic current vanishes 
[23] . Similarly, as discussed in Ref. [13], there is no leading renormalon contribution to this 
operator, which implies + = in the renormalon model. In App. [C] we carefully 
investigate a number of different sum rules for ± which are mutually consistent and 
yield the results given in Tab. [2], with C^ + C^- ^ 0. On the other hand, our result for — 
is consistent with older estimates based on the analysis of the leading contribution to the 
chiral-odd sum rules [23j, albeit a factor two smaller. 

The resulting four two-particle twist-4 DAs, as given by ([S2SD, dSSSj), and (H^Sj), 

are shown in Figs. [1] and [2l We use the renormalon- model predictions for all matrix elements 
which are not known from a direct calculation. For C4^'^, we use ^4^'^ = (^4^ — Ct)/'^- The 
SU(3) breaking is moderate in ^f"*", but obviously more pronounced for ip\^ . This feature is 
mainly due to the different shape of the asymptotic DAs which vanish at the end-points for 
04, but are non-zero for ?/'4. As is seen from the behaviour of iii particular. Fig. [T], the 
finite mass corrections in ms change the shape of the DA noticeably for u — > 1; this result is 
dominated by the terms linear in m^. Keeping all quark masses, the behaviour very close to 
the end-points is given by mq{ms — mq) Inu for m — > 1 and ms{ms — rriq) Inu for m — >■ 0. For 
M — i> the logarithmic term is dominant and causes the marked asymmetry in the dashed 
(green) curve in the left panel of Fig. [H For the (j) meson, the logarithms vanish as rriq —>■ rris- 
A similar effect is at play for , Fig. [21 but is slightly less marked numerically. Due to 
the dominance of these finite-mass corrections, the renormalon model alone gives only a 
poor description of the full DAs, see Fig. [31 In particular the size of the apex at m = 1/2 is 
considerably underestimated. The apex is actually dominated by the contribution of ajj' ^ = 1 
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Figure 1: [Colour online] Left panel: ipl, ()3.26p . as a function of u for the central value of the 
hadronic parameters, for fj, = 1 GeV. Solid [red] line: ipl.p, dashed [green]: ip^.j^,, short-dashed 
[blue[: tp\.^- Right panel: same for (/)|, (j3.28p . 
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Figure 2: [Colour onhne] Same as Fig. [D for V'f, (Oil) , and 0f , (05]) . 
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Figure 3: [Colour online] Renormalon-model predictions for 0^;^ and V^lpi '^i-p = in the renor- 
malon model. 
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to the DAs, see Eqs. (13.271) and (13.291) . Technically speaking, this term is a mass correction 
and hence not included in the renormalon model. Despite this shortcoming, the renormalon 
model is very useful for estimating otherwise only poorly constrained higher- conf or mal waves 
of twist-4 DAs, and in particular G-parity-breaking parameters. 

6 Summary and Conclusions 

In this paper we have studied the twist-4 two- and three-particle distribution amplitudes 
of p, K* and mesons in QCD and expressed them in a model-independent way by a 
minimal number of non-perturbative parameters. The work presented here is an extension 
of Refs. [151 [161 [17] ci-nd completes the analysis of SU(3)-breaking corrections by also including 
G-parity-breaking corrections in irts—niq to twist-4 distribution amplitudes. Our main results 
are the expressions for twist-4 two-particle distribution amphtudes, Eqs. (I3.26p . (I3.28p . (14.241) . 
(14.251) . and the complete set of twist-4 input parameters. Tab. [2l With these results, a 
complete set of light-meson DAs of twist 2, 3 and 4 is available for both pseudoscalar and 
vector mesons. 

Our approach consists of two components. One is the use of the QCD equations of 
motion, which allow dynamically dependent DAs to be expressed in terms of independent 
ones. The other ingredient is conformal expansion, which makes it possible to separate 
transverse and longitudinal variables in the wave functions, the former ones being governed 
by renormalisation-group equations, the latter ones being described in terms of irreducible 
representations of the corresponding symmetry group. We have derived expressions for all 
twist-4 two- and three-particle distribution amplitudes to NLO in the conformal expansion, 
including both chiral corrections 0{ms + mq) and G-parity-breaking corrections 0{ms — mq)] 
the corresponding formulas are given in Sees. [3] and HI We have also generalized the renor- 
malon model of Ref. [13] to describe SU(3)-breaking contributions to high-order conformal 
partial waves. 

We have done a complete reanalysis of the numerical values of the relevant higher-twist 
hadronic parameters from QCD sum rules. Our sum rules can be compared, in the chiral 
limit, with existing calculations for the p [161 [2Z|- Whenever possible, we have aimed at 
determining these matrix elements from more than one sum rule; we find mutually consistent 
results, which provides a consistency check of the approach. Our final numerical results, at 
the scales 1 and 2 GeV, are collected in Tab. [2l Any substantial improvement of these 
results will require input from alternative non-perturbative methods, in particular lattice 
calculations. 

We hope that our results will contribute to a better understanding of SU(3)-breaking 
effects in hard exclusive processes and in particular in the decays of B^^d and Bg mesons into 
final states containing light vector mesons. 
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Appendices 

A Non-Local Operator Identities 

For completeness, we quote the following non-local operator identities from Refs. [TSl EI] : 

d _ 

— q{x)'-i^s{-x) = J dvvq{x)xagGa^,{vx)'yf,s{-x) 

+i{ms + mg)q{x)s{—x), (A.l) 
d^{q{x)'j^s{-x)} = -i J dvq{x)xagGaf,{vx)'j^s{-x) 

- i{ms- mq)q{x)s{-x), (A.2) 

d _ - 
d^,q{x)(y^^,s{-x) = -i — q{x)s{-x) + J dvvq{x)xpgG py{vx)s{-x) 

-i j dvq{x)xpgGpp{vx)apys{-x) 

-{rris + mq)g(x)7^s(-x), (A.3) 

d 

q{x)(Tpus{—x) = —iduq{x)s{—x) + / dv q{x)xpgG pu{vx)s{—x) 



dxp j_i 



— 2 y dvvq{x)XpgGpp{vx)o^ys{—x) 

+ {rris - mg)q{x)'y^s{-x). (A.4) 



Here dp is the total derivative defined as 



d_ 



df, {q{x)Ts{-x)} = t;— {q{x + y)[x + y, -x + y]Ts{-x + y)} 



By taking matrix elements of the above relations between the vacuum and the meson state, 
one obtains exact integral representations for those DAs that are not dynamically indepen- 
dent. 
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B Chiral-Even Twist-4 Parameters 



In this appendix we calculate the parameters Qy and ujly defined in Eqs. (13.121) and (13.191) . 
We shall obtain them from QCD sum rules, using various correlation functions of either 
identical currents, so-called diagonal sum rules, or different currents, so-called non-diagonal 
sum rules. We shall further distinguish between pure-parity (PP) and mixed-parity (MP) 
sum rules, depending on the parity of hadronic states that contribute to these correlation 
functions. 

Let us first consider the non-diagonal correlation function 

z^z'^t J dSe-'^^Q\Tq{z)gG,^{vz)^''^,s{my)luq{ym = {vzf j Vae-'^P^^"^+''-'\{a), 

(B.l) 

where both currents have the same chirality. The integration measure Va is defined in (12.171) . 
This is a MP correlation function, with both = 1~ and 0"'' states contributing. We have 
calculated the OPE including condensates up to dimension 6: 



1 



+ 



1 



—p In — ^aia2a3. , . , , . , 2 \ 
zir-^ —p^ (^ai(l — tti) ^2(1 — ^2) J W \ 



1 /a 



1 Ots 
3p^ TT 

1 as 



mq{qq) as^as - 3)6{a2) - a3a3(5'(a2) 



3p^ TT 

2 a. 



+ 



3p^ TT 
1 



/i 



^3 In - 0^3 ln(Q;3a3) + 0:303 



-p- 



p4 



1 



—TTas{ss)^ + ^ms{sagGs) 5(01)5(03) 
8 1 

—TTas{qq)^ + -mq{qagGq) 5(02)^(03) 



(B.2) 



16 

- — 7ro,(gg)(ss)5(oi)(5(o2) . 

In the local limit and zero quark masses, the result agrees with the calculation in Ref. 
In the product of 5-functions 5'(oi_2)5(oi + 02 + 03 — 1), 6 has to be integrated over before 
6'. 

From flB.2p we obtain the following sum rules: 



--^ M'{1- r(2, .o/M^) ] + \- i^M) + ms{-ss)) 
+ 1 — {rusiqq) + mg{ss)) + 7^; - In ^ + 1(0, Sq/M' 
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(-0.24 ± 0.01)^ GeV^ 


(ss) = (1 - 53) (qq) 


{qagGq) = 


m^ {qq) 


{sagGs) = (1 - d^JiqagGq) 




(0.012 ± 0.006) GeV^ 


(gy-G^) = (0.08 ± 0.02) GeV^ [28j 


ml 


= (0.8±0.1)GeV^ 6 


3 = 0.2 ±0.2, ^5 = 0.2 ±0.2 



m,(2 GeV) = (100 ± 20) MeV 



m,(l GeV) = (133 ± 27) MeV 



a^(m^) = 0.1176 ±0.002 < — > a,(l GeV) = 0.497 ± 0.005 

Table A: Input parameters for sum rules at the renormalization scale /U = 1 GeV. The value of nig 
is obtained from unquenched lattice calculations with Nf = 2 light quark flavours as summarized in 
Ref. [29], which agrees with the results from QCD sum-rule calculations [30]. as{mz) is the PDG 
average [3T] . 



6 \ vr 



3M2 



[mq{qgGaq) + ms{sgGas)) 



27M2 

Sets 
2592^ 

11 a 



9M2 



(B.3) 



19 (T 

M'{i- r(2, VM')} - 7^ ^ K(^g) ± nis{-ss)) 



648 vr 



324 TT 



M2 



ms{qq) + mg{ss)) <^ - ± 7^ - In — ± 1(0, Sq/M 



where 
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27M2 



{mq{qgGaq) + ms{sgGas)) 



{{qq?+{ssf]-TTT7^m{ss) 



243M2 
r(a,So/M2) 



81M2 



(B.4) 



(1B.3P follows from (1B.2|) by integration over Pa with the weight factor 1, and (IB. 40 by 
integration with weight factor (03 — 4/9)/2, see Eq. f l3.19p . The sum rules for p are obtained 
by letting s q and those for (p by letting q s. 

We evaluate the above sum rules using the input parameters collected in Tabs. [1] and \^ 
the results, for central values of the input parameters, are shown in Fig. \^ The sum rules 
are dominated by the contribution of the gluon condensate; the impact of the specific value 
of the continuum threshold sq is only moderate. The figure also shows that the impact of 
SU(3)-breaking is very relevant: the contributions of the quark and mixed condensate reduce 
the values of Cj^ ^ind uj\y In the Bor el- window 1 GeV^ < < 2 GeV^, and including the 
input-parameter uncertainties given in Tab. [Aj we find, at the scale /x = 1 GeV: 



4p 



0.07 ±0.03 



UJ 



4p 



-0.03 ±0.01 
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[GeV'] [GeV^] 

Figure A: [Colour online] cly (left) and uj\y (right) from the non-diagonal MP sum rules ()B.3P 
and ()B.4p as functions of M'^, for central values of the input parameters. Solid [red] lines: p 
{so = 1.5 GeV^), long dashes [green]: K* {sq = 1.8 GeV^), short dashes [blue]: (p {sq = 2GeV'^). 
All parameters are evaluated at the scale /i = 1 GeV. 



= 0.02 ± 0.02 , = -0.02 ± 0.01 



Q = 0.00 ±0.02, =-0.02 ±0.01. (B.5) 



We have added in quadrature all individual sources of uncertainty. The total error is domi- 
nated by that of the gluon condensate. 

For (li^, , we also consider diagonal sum rules which can be obtained from the correlation 
function 

K = ^ jd'xe^^^{^\TJl{x){j'^)\m = P.P.Kip") - 9,.T^X{p') , (B.6) 

with the current = q gG ^alal^s. For p, 11^^ was calculated in Ref. |2^, while the 
SU(3)-corrections were calculated in Ref. [H], including contributions from condensates up to 
dimension 8. The suitability of this correlation function for extracting (^|^* is not immediately 
obvious: 11^ contains contributions not only of vector mesons, but also of hybrid O"*" mesons, 
whose coupling to is much larger than that of the K*, ruling out the possibility to 
construct a MP sum rule for C\k*- This situation is qualitatively different from that of the 
non-diagonal correlation function fIB.ll) . where the presence of s'^yQ removes all contributions 
from hybrid mesons. We hence focus on the PP function II]^. From Ref. we quote 



2407r3 ^ -p2 36 \ TT / -p 



,2 



2 2 
r /-\l2i 1^ I r /-\, /-\l2i A'' 



+ —[m„{qq) + ms{ss)]p In 5- ± 77— [m^(gg) ± mg(ss)]p In 

DTT — p"^ ISvr — 



9 ^^^'^ ' 1927r2 ' 
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[f^q{<l<^gGq) + ms{sagGs)] In ■ ^ 



^44 ^ L ^ ^/ ■ ^\ ^ /J 

^ [ms{qagGq) + mq{sagGs)] In ■ ^ 



144 ^ ''^^ ^ ^' ' ^ 'J -p2 

The PP sum rule for (^|^* is 

(/i)'<.(CLO'e-™-/*^' =Ssubnr, (B.8) 
where Sgubn]^ is the continuum-subtracted Borel transform of II^ , which we define as 

Jo S — p Jq 

in terms of the dispersive representation of HY ■ For p, the above sum rule was derived and 
analysed in Ref. [27]. It features a large negative contribution from the gluon condensate 
which, for Borel parameters ~ 2 GeV^ and continuum thresholds sq between 1.3 and 
3 GeV^, drives the right-hand side of fIB.Sp negative. In Ref. it was argued that this 



large negative contribution signals the presence of a larger mass scale ~ 2 GeV in the spectral 
density and is indicative for a breakdown of quark-hadron duality, at least if the usual simple 
continuum model with only one resonance, the p, is used. A remedy is the use of a more 
appropriate continuum model including higher mass states like p(1450). This automatically 
increases Sq, but it also turns out that the coupling of p(1450) to the gluonic current is 
larger than that of p(770), which does not really help the determination of (Ip. An alternative 
is to analyse (IB.SP for small ^ IGeV^, where duality still works and the suppression 
of higher-mass resonances is effective. Numerically, the sum rule is then dominated by the 
gluon and the dimension-8 condensate {qq){qcrgGq). Clearly such a sum rule cannot give an 
accurate estimate of (ly, so we only use it as a consistency check for the results obtained 



from (IB.3P and, in particular, the large SU(3) breaking. The results from fIB.SP are shown 
in Fig. [HI Note that the breakdown of duality sets in the earlier, the heavier the meson. 
Although it is not possible to extract precise values for (ly, we see that the values are not 
inconsistent with the results from the non-diagonal sum rule, f IB.Sp . and that in particular 
the relative hierarchy, (f^ > ^|^, > is reproduced. 

Our final results are given in (IB.Sp . A comparison with earlier determinations is given in 
Sec. El 
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Figure B: [Colour online] Qy from (|B.8p . Solid [red] line: p, long dashes [green]: K* , short dashes 
[blue]: (j). Same input parameters as in Fig. lAl 

C Chiral-odd Twist-4 Parameters 

In this appendix we calculate 

Ci = C4V±ci.*- (c.i) 

Like for C,\k*^ consider both non-diagonal and diagonal sum rules - the former for all 
mesons, the latter only for p. We also consider PP and MP sum rules. To distinguish 
between the results of these sum rules, the following notation proves convenient: 



^± Id(nd),pp(mp) ■ 

Let us start with the non-diagonal sum rules for (^^, yielding C,^ 
correlation function is 



ND,PP(MP) 



(C.2) 
. The relevant 



= ^ j Ae-^^^(0|r[g(7(G^.±^G^.75)5](0)[s-a„^g](y)|0) 
= ^{ [{Pfi9uaPp) - (/^ ^ z^)] - [a ^ /?] 

The invariant functions Hy and contain contributions of 1~ and 1^ mesons, respectively, 
and can be separated by considering the two projections 



n± + n± 



17='= = n^^'^n^l^W^ — — 



In calculating 11^2^ '^^ dimensional regularization and the identity 



(C.4) 



(C.5) 
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in order to avoid ambiguities with the definition of the epsilon-tensor and the 75 matrix in d 
dimensions. Here {...,...} denotes the anti-commutator. We find 



n 



V 



n 



V 



487r3 
1 



In 



1 a. 



-p^ 3 vr 



+ 



I2p 

as 



ms{qq) + mq{ss) 
1 



In 



IT 

^p2 Qp2 



mq{qagGq) + ms{sagGs) 



+ 



mq{sagGs) + nigiqcrgGq) 



a. 



jns{qq) + mq{ss)] - — [ms{ss) + mg{qq)] 
yTT ovr 



a. 



1 /^^ la. 
ryw m H 



+ 



12p2 

a 



TUsiqq) + mg{ss) 
1 



mq{qagGq) + ms{sagGs) 



6p2 



_pZ gpZ 24 \ TT 



mq{sagGs) + ms{qagGq) 



a. 



,ms{qq) + mq{ss)] - [ms{ss) + mq{qq)] 

yTT DVT 



727r3 



j9 In 



— p2 9 TT 



12 \ TT 



In 



+ 



1 «^ 
9 vr 

1 



In- 



mq{qagGq) + ms{sagGs) 



p2 gp2 

5a 



— 

\2 I /;^„\2 



,2 

■ 12 \ TT 



.__p in_^ + _(_G )ln- 



27n 



27p2 

rrisiqq) + mg{ss)] 



{ssy + {qqY 



—p^ 9 TT 



m^(gg) + mq{ss) 



In 



1 a, 
9^ L 

1 

TV - 



mq{qq) + ms(ss) 



In ^ - -—ras{ss){qq) - ^r^a 



mq{qagGq) + ms{sagGs) 



p2 Qp2 

5a« 



277r 



27p2 



—p^ 
{ss)^ + {qqf 



(C.6) 



Again, the correlation functions for the p meson are obtained by s g, and those for by 
q ^ s. The above functions allow one to derive the PP sum rules 



3±InD,PP'^ " — '^sub J-J-v , (C.7) 

and correspondingly for axial- vector mesons. As discussed in Ref. [20], there are actually 
two strange l"*" mesons, Ki{127Q) and i^'i(1400), which are usually interpreted as mixture of 
a ^Pi state, the Ka-, and a ^Pi state, the K^, [32[ [33] : 

i\:i(1270) = KaCOsdK - KbSiYlOK, 

Ki{imQ) = KaSineK + KbCosOK. 
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Figure C: [Colour online] C:j^|j^j^pp (left) and C-|n£,pp (right) from the PP sum rule (|C.7p . as 
functions of M^, for fj. = 1 GeV and central values of input parameters. Solid [red] lines: p 
{so = 1.2 GeV^), long dashes [green]: K* (sq = 1.4 GeV^), short dashes [blue]: (/> (sq = 1.8 GeV^). 

The results of Refs. [32l[33] indicate that the system is close to ideal mixing, i.e. 9k ~ 45°. To 
the accuracy needed here it is then sufficient to replace the two resonances by one effective one 
with the mass ttiki = 1.34 GeV [33j. We find that the sum rules for and its axial- vector 
equivalent are dominated by the gluon condensate contribution, which implies 

C+(i^i)|ND,PP ^ C+Ld,pp (C.8) 

with ~ 30% accuracy. For no single contribution is dominant, but one still finds that 
and (^{Ki) have opposite sign. This is similar to the situation with PP sum rules for 
the G-odd twist-4 parameters /tf^. and k;^^^ discussed in Ref . [2IJ , and has some impact on 
the relative size of continuum contributions in PP vs. MP sum rules. We will come back to 
this point below. In Fig. \C\ we show the results for Cil^opp- The values of sq are chosen 
in such a way as to ensure maximum stability of ^jj'rjj^j^pp in the Borel parameter. Note 
that the sum rules (1G.7P are quite sensitive to the value of the continuum threshold, which 
we vary by ±0.3 GeV^. Including this uncertainty and the variation in M^, in the interval 
1 GeV^ < M2 < 2 GeV^ and the error induced by the hadronic input parameters. Tabs. |2] 
and El we find, at the scale = 1 GeV: 

C+(p)|nd,pp = -0-12 ±0.04, C-(p)Ld,pp = 0.03 ± 0.02 , 
C^(ir*)|^j3 PP = -0.07 ±0.02, C-(^*)Ld,pp = 0-04 ±0.02, 

C^(0)|^j3PP = -0.05 ±0.02, C-(0)Ld,pp =0.04±0.02. (C.9) 

We have added all individual uncertainties in quadrature. The bulk of SU(3) breaking in 
these quantities is due to the factor my{fyY in (IC.7I) . with m^^lf^,)"^ / {mp{f^)'^) = 2.2 
and rn'^{f^y/{mp{fp)'^) = 3.8, which explains the relative hierarchy |C+(p)| > |C+(-^*)l > 
|C:jl:(0)|. For one has a cancellation of several terms which renders the interpretation of 
the hierarchy of the curves in Fig. Oless clear-cut. 
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Figure D: [Colour online] C+|ndmp 0^^^^) ^^"^ ^-Indmp ("g^*) ^^'^ sum rule (jC.lOp . 

as functions of M^, for /i = 1 GeV and central values of input parameters. Solid [red] lines: p 
{so = 1.0 GeV^), long dashes [green]: K* (sq = 1.2 GeV^), short dashes [blue]: cf) {sq = 1.5 GeV^). 



Let us now turn to the extraction of C± from MP sum rules, which contain contributions 
from both vector and axial- vector mesons. These sum rules are derived from the combinations 
(Ily + n^)/p^ and, thanks to the factor benefit from a smaller mass dimension than 

the PP sum rules: 



m. 



K* S± |nD,MP 



1 



i3,,b- (n^ + n±). 



(C.IO) 



consists of 



The corresponding results are shown in Fig. [Dl The sum rule for C+In^mp 
only two terms: the quark-condensate and the four-quark-condensate contributions, with 
the latter dominant. Such a sum rule, sensitive to higher-dimensional condensates, is not 
reliable and we do not include its results into our final value for C- Ij^^ y^^p, on the other 
hand, does not receive any contribution from the four-quark condensate, and is dominated by 
the gluon condensate. In contrast to Cil^opp; C±|ndmp rather insensitive to the precise 
value of the continuum threshold Sq. The reason for this is the different sign of vector and 
axial-vector contributions to the hadronic side of the sum rule: as mentioned above, the 1+ 
and 1~ matrix elements tend to have different sign, and hence the resonance contributions 
to the sum rule tend to cancel, reducing the size of the continuum contribution. We choose 
So for the MP sum rules slightly below that for PP sum rules, to account for the lower mass 
of the 1^ ground state as compared to the first 1~ excitation. Using again the hadronic 
input parameters from Tab. [2] and Tab. \^ and varying sq by ±0.3 GeV^ and in the 
window 1 GeV^ < < 2 GeV^, we obtain the following results for C- |ndmp scale 
/i = lGeV): 



(Hp) 



ND.MP 



0.07 ±0.03, C-(^*) 



IND,MP 



0.03 ±0.02, C^(0) 



ND,MP 



0±0.02. 
(C.ll) 

We have added all individual uncertainties in quadrature. We do not give any results for 
C+ I p^j^ because this sum rule is dominated by higher-dimension condensates and hence 
not reliable. Again SU(3) breaking in the above numbers is dominated by the overall factors 
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m 



(/y)^ in (]C.10|) . Both PP and MP non-diagonal sum rules agree about the signs of C±: 



C > and C+ < 0. 

Let us now discuss the diagonal sum rules for which can be derived from the correlation 
function of two quark- ant iquark-gluon currents: 

^aiu = ^ [ Ae-^^^(0|T{[gr7(G^,±^G^,75)5](0)[^^?(G,;3±zG,^75)g](2/)}|0) 



Like for the non-diagonal correlation function (IC.Sp . the invariant functions and 
contain contributions of = 1~ and 1"^ states, respectively, and can be separated by the 
projections 

^{pzf 



2tT±± 



g ^^aBh 



n±± ^ g'-g^'n^p,^ = -6 



(C.13) 



Obviously H^^^^-j = H^^-j. We calculate these invariant functions neglecting mass corrections 
and find 
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(C.14) 
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Figure E: [Colour online] (C+C-)/) ^ov p from the PP sum rule ()C.16p (solid [red] curve) and the 
MP sum rule (IC.lSj) (dashed [green] curve). 
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Figure F: [Colour online] |C+|d (left) and |C-|d (right), for the p, from the sum rules (|C.15p 
and ()C.17p . Solid [red] curves: PP (sq = 1.2 GeV^), dashed [green] curves: MP sum rules (sq = 
1.0 GeV^). 



Without quark-mass corrections, the above results only allow the calculation of the p cou- 
plings. We construct PP and MP sum rules for (C+)^ and (C-)^, and also for the product 



PP 



D,MP 



f2 77='=='= 
*-'sub-'--'-y 5 



i3,,b- (n±± + n±±) 



p 



(C+e)|DMP = ^-b-^ (n^ + nr) . 



(C.15) 
(C.16) 

(C.17) 
(C.18) 



The results are shown in Figs. [E] and [0 As Fig. [E] indicates, both flC.161) and flC.lSP predict 
different sign for and C^, which agrees with the result from non-diagonal sum rules, Figs.lCl 
and [Dl In Fig. [F] we plot the results from (IC.ISP and ( 1C.17I) . All four sum rules receive 
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contributions only from perturbation theory and the dimension 8 condensate {qq){q(TgGq). 
For both sum rules for perturbation theory is dominant, whereas for those for the 
condensate is dominant. For C-|dmp comes with a negative sign, which explains the fact 
that the dashed curve in Fig. [F] (right) only starts at large M^: below that, the result is 
imaginary. A possible explanation is the contribution of hybrid 1^ states to the sum rules, 
both for C- lo jy[p and C- pp- Indeed, replacing the parametrisation on the left-hand side of 
flC.171) by a contribution of such a state, with a mass ~ 1.5 or 2 GeV, and correspondingly 
larger sq ~ 4 GeV^, the coupling becomes real. In any case, the dominance of the dimension- 
8 term in (qq) (qagGq) in these sum rules renders them unreliable and hence we discard their 
results. As for C+ Ip,; the sum rule is very dependent on sq, which totally dominates the error 
budget. This is another manifestation of the fact that the currents in (]C.12I) have a strong 
coupling to hybrid states. We finally find 

C+ Id,pp = -0-06 ± 0.02, C+ |d,mp = '^-^^ ± 0-02 • (C.19) 

Our last task is to give final results for and, equivalently, Civ and according to 
(]C.1|) . Let us first discuss the p parameters, as only for those we have information from both 



diagonal and non-diagonal sum rules. For we have three results, C+|p,pp and C+lpi^p 
from flC.191) . dominated by perturbation theory, and ^jj'rjj^j^pp = —0.12 ± 0.04, dominated 
by the gluon condensate contribution. As none of these results is a priori "better" than the 
others, we average over all of them to obtain our final result 

C+(p, 1 GeV) = -0.10 ±0.06. (C.20) 

The average is smaller than the MP sum result alone, which we take into account when 
arriving at our final results for K* and from (1G.9P : 

C+ (i^*,l GeV) = -0.06 ±0.03, C+ (0, 1 GeV) = -0.04 ± 0.03 . (C.21) 

As for the diagonal sum rules have to be discarded, whereas the non-diagonal ones yield 
^- Indpp ~ ^'^"^ 0-02, with the most relevant contributions from perturbation theory and 
the dimension 6 condensate (gg)^, and C-I^pi^p ~ 0.07 ± 0.03, with the most relevant 
contributions from perturbation theory and the gluon condensate. Here we obtain our final 
result as a straight average of PP and MP sum rules and find, from ( ]C.9P and (IC.lip : 



Cf(p,lGeV) = 0.05±0.05, C-(^*,lGeV) = 0.04±0.04, C-(0,lGeV) = 0.02±0.04. 

(C.22) 

From (IC.ip . we also find the final results for (^y and Civ- 



Ci(l GeV) = -0.03 ± 0.05 , (4^(1 GeV) = -0.01 ± 0.03 , C4^^(l GeV) = -0.01 ± 0.03 , 

Ci(lGeV) = -0.08 ±0.05, Ci^.(l GeV) = -0.05 ± 0.04 , Ci(l GeV) = -0.03 ± 0.04 . 

(C.23) 

A comparison of our results with those from previous calculations is given in Sec. O 
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